Optical properties and electronic structure of a'-Nai-zCa^A^Os 



O 

o 
o 



C. Presura 1 , D. van der Marel 1 , M. Dischner 2 , C. Geibel 2 , R.K. Kremer 3 
1 Solid State Physics Laboratory, University of Groningen, Nijenborgh 4, 9747 AG Groningen, The Netherlands 
2 Max- Planck- Institut fur Chemische Physik fester Stoffe, D-01187 Dresden, Germany 
3 Max- Planck- Institut fur Festkorperforschung, Heisenbergstrasse 1, 70569 Stuttgart, Germany 

(February 1, 2008) 

The dielectric function of a'-Nai-zCa^ V2O5 (0 < x < 20%) was measured for the a and b axes 
in the photon energy range 0.8-4.5 eV at room temperature. By varying the Ca-concentration we 
control the relative abundancy of V 4+ and V 5+ . We observe that the intensity of the main optical 
absorption peak at 1 eV is proportional to the number of V 5+ -ions. This rules out the interpretation 
as a V 4+ d — d excitation, and it establishes that this is the on-rung bonding-antibonding transition. 

PACS numbers: 78.40.-q, 71.35.-y, 75.50.-y 
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I. INTRODUCTION 

a'-Nai_ a; Ca a: V205 belongs to the larger group of a' 
vanadium pentoxides, with-the chemical formula AV2O5 
(A= Li, Na, Ca, Mg, etc.)tl. Their structure is remark- 
ably similar to that of the parent compound V2O5, which 
consists of layers of square pyramids of O surrounding a 
V +5 ion. The layers are kept together via weak forces, 
which account for the easy cleavage of this oxide along 
(001). The basic building blocks forming the V2O5 lay- 
ers are parallel ladders of VO5 pyramids. In AV2O5 the 
A atoms enter the space between the layers and act as 
electron donors for the V2O5 layers. In the case of a'- 
NaV2C>5, every doped electron is shared between two V 
atoms. As a result the average valence of the V-ions 
corresponds to V +4 5 . X-ray diffraction indicates that 
at room temperature all V-ions are in the same mixed 
valence statdoQ. 

Partial substitution of the Na + with Ca 2+ leaves the 
a' crystal structure intact, but alters the relative abun- 
dance of V 4 + and V 5+ , N 4 + : N 5+ = (1 + x) : (1 - x). 
In this paper we report spectroscopic ellipsometry mea- 
surements on a'-Nai_ a: Ca a ;V205 (x=0, 0.06, 0.15 and 
0.20), in the energy range 0.8-4.5 eV. We employ the de- 
pendence of the optical spectra on the A^ 4+ : iV 5+ ratio 
to identify the main components in the optical spectra, 
which in turn we use to reveal the electronic structure of 
this material. 



II. DETAILS OF SAMPLE PREPARATION AND 
EXPERIMENTAL SETUP 

The crystals (CR8, 45008, and 45010) had dimensions 
of approximately 2, 1 and 0.3 mm along the a, b, and 
c axes respectively. The samples 45008 and 45010 were 
prepared from NaV03 fluxp. In a first step a mixture of 
Na 2 C03 and V2O5 is heated up to 550 C in air to form 
NaVC>3. In a second step the NaVC>3 is mixed with VO2 
in the ratio of 8:1 and then heated up to 800 C in an evac- 
uated quartz tube and cooled down at a rate of 1 K per 



hour. The excess NaVC>3 was dissolved in water. Then 
the doped samples were produced by substituting in the 
first step Na2CC>3 by CaC03. The chemical composi- 
tion of the samples has been determined using Energy 
Dispersive X-ray Fluoresence microprobe measurements. 
The results showed that the real Ca content of some sam- 
ples was smaller that the nominal one (with a factor of 
0.75), and that position dependent variations of the Na 
stoichiometry are below 2%. A standard spectroscopic 
ellipsometer was used to collect ellipsometric data in the 
range of 6000 to 35000 cm— 1 from the ab planes of the 
crystals using two different crystal orientations, and to 
measure normal incidence reflectivity spectra of the ac 
plane with the electric field vector along the c-direction. 



III. DATA COLLECTION AND ANALYZES 

We performed ellipsometric measurements on the (001) 
surfaces of the crystals both with the plane of incidence 
of the light along the a and the b axis. An angle of inci- 
dence 9, of 66°, was used in all experiments. Ellipsometry 
provides directly the amplitude and phase of the ratio 
the reflectivity coefficients of p- and s-polarized lighi 
r p (ui)/r s (ijj). For an anisotropic crystal with the three 
optical axes ariented along the surface normal (p _L), 
perpendicular to the plane of incidence (s), and along 
the intersection of the plane of incidence and the surface 
(p ||), this ratio is related to the dielectric tensor elements 



along these three directions (e p j_, 
to the expression: 
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To extract the dielectric constant from the ellipsometric 
parameters we proceed in two steps: First the pseudo- 
dielectric functions along the optical axes are extracted 
from the ellipsometric data using the inversion formula 
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FIG. 1. (a) Pseude-optical conductivity a\ s {uj) of 
«'-NaV 2 5 at T = 300KB. The spectra were taken on the 
(001) surface with a and b axes successively lying in the plane 
of incidence, (b) E\\c reflectivity at T = 300K. (c) Optical 
conductivity a a {uj) and <Jb{u>) corrected for contributions of 
the c-axis dielectric function to af a (ui). 



For isotropic crystals this expression provides the dielec- 
tric function directly. The pseudo-dielectric function is 
close to the dielectric tensor elements along the intersea- 
tion of the plane of incidence and the crystal surface!! 
A biaxial crystal like a'-NaV20s has three complex di- 
electric functions, e a , e& and e c along each optical axis, 
and an ellipsometric measurement involves all three ten- 
sor components of the dielectric matrix. In addition to 
the pseudo-dielectric functions displayed Fig.[|a, e c (oj) is 
required. No ac or be crystal-planes were available large 
enough to do ellipsometry with our setup. We there- 
fore measured the c-axis reflectivity (Fig. Hp) of the bc- 
plane of the pristine material (sample CR3). The spec- 
trum contains no (or very weak) absorption peaks in-,thc 
measured frequency range, as it was reported earlierBE2l, 
providing a very reliable determination of the dielectric 
function e c using Kramers-Kronig analysis. Due to the 
absence of strong resonances, e c has a minor influence on 
the recorded ellipsometric spectra. In Fig. 0c the optical 
conductivity is displayed taking into account all correc- 
tions due to the anisotropy. We see that the conversion 
from e ps (u>) to e(ui) has indeed a rather small effect on 
the spectra. In essence it leads to a factor 0.5 re-scaling 
of ct(uj). 

The are in general agreement with previous 

resultsfcyilJ using Kramers-Kroning analysis of reflectiv- 
ity data. Along the a-direction we observe a peak at 0.9 
eV with a shoulder at 1.4 eV, a peak at 3.3 eV and the 
slope of a peak above 4.2 eV, outside our spectral win- 
dow. A similar blue-shifted sequence is observed along 
the b-direction. The leV peak drops rather sharply and 
extrapolates to zero at 0.7 eV. However, weak absorption 
has been observed within the entire far and mid-infrared 
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FIG. 2. Optical conductivity at T = 300K of Re cr(w) of 
Ca^Nai-^Os for x=0 (solid), x=0.15 (dash) and x=0.20 
(dot): E || a (panel a) and E \\ b (panel b) 
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FIG. 3. Intensity of peak A calculated by integrating opti- 
cal conductivity up to 1.9eV (squares), or by using the height 
of the peak (circles) in Ca^Nai-zX^Os for x=0.06, x=0.15 
and x=20 and in Na3;V205 for x=0.95. The dashed line indi- 
cates the theoretical intensity of in-site Vriri transitions versus 
doping. The dotted line indicates the theoretical intensity ac- 
cording to the model Hamiltonian Eq.H. 
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rang The strong optical absorbtion within the en- 

tire visible spectrum causes the characteristic black ap- 
pearance of this material. 

The peak positions appear to be doping independent, but 
the striking observation is that the intensity of the peaks 
depends strongly on doping. In particular, the measured 
intensity of the 1 eV peak for the a axis is directly pro- 
portional to 1 — x (Fig. U). The 1 eV peak for the b axis 
shows also a decrease upon doping. 



IV. MAIN ELEMENTS OF THE ELECTRONIC 
STRUCTURE 

Before entering the interpretation of the data, we need 
to discuss in some more detail the main elements of the 
electronic structure of these compounds. The basic build- 
ing block of the crystal structure of a'-NaA^Os is formed 
by VOV dimers. These dimers form the rungs of quasi 
one-dimensional ladders. The V-ions forming the rungs 
are bonded along the ladder direction via oxygen ions. 
The backbone of the electronic structure is formed by the 
oxygen 2p and V3<i states. Photoelectron spectroscopjlij 
has provided crucial information on the occupied elec- 
tronic levels: The oxygen 2p states have the lowest en- 
ergy (the highest binding energy). They form a band 
about 4 eV wide, which is fully occupied. The occupied 
part of the V 3d states is located about 3 eV above the 
top of the oxygen bands. Due to ligand field splittings 
the V 3d manifold is spread over a range of at least 3 
eV. The 3d xy is of the V 4+ is occupied with one elec- 
tron. The unoccupied d xz and d yz levels have an energy 
at least 1 eV higher. These in turn areJocated about 2 
eV below the d x 2_ y 2 and the d z i levelscl. The relevance 
of the O 2p bands is that they provide a path for virtual 
hopping processes between the V-sites. The coupling be- 
tween V-sites is through virtual hopping via 7r-bonded O 
2p y states, indicated schematically in FigJ^D. The effec- 
tive hopping parameter between V-sites is t± — 0.3eV 
on tha-saiae-rung, and t± » — 0.2eV along the legs of the 
laddeni-aOO. The number of electrons is one per pair of 
V-atoms. Approaching the ladders as a linear array of 
rungs, weakly coupled along the direction of the ladder, 
results in a model of electrons occupying a narrow band 
of states formed by the anti-symmetric combination of 
the two V 3d-states forming the rungs, hereafter referred 
to as V-V bonding levels. 

Hence, we see that the basic building block are the pairs 
of V3d states, together forming the rungs of the ladders. 
The essential charge and spin degrees of freedom of a sin- 
gle rung are identical to the Heitler-London model of the 
H^ ion, pith the V3d xy states playing the role of the H 
Is statestil. The relevant Hamiltonion is 
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FIG. 4. Diagram of the many-body eigenstates of a VOV 
rung, occupied with 0, 1 and 2 electrons, obtained within 
the Heitler-London model for U = 4eV A = 0.8eV, and and 
t± = 0.3eV. 



where d' L , R *. a crea tes an electron in the lefthand (right- 
hand) d xy orbital on the rung, and a is the spin-index. 
The bias potential A between the two V-sites accounts for 
a possible left/right charge imbalance!! The most rele- 
vant states for the ground state are d^ xy and dji xy . Pure 
NaV2C>5 contains one electron per rung in the ground 
state. Important in the present discussion are the eigen- 
states and energies of a rung with 0, 1 or 2 electrons. The 
eigenstates and energies are listed in Table |. 
In Fig. [I| the level diagram is displayed. In this rep- 
resentation N e — 0(2) corresponds to the one electron 
removal(addition) states, for noninteracting electron pic- 
ture indicated as the "occupied" ("empty") states. In 
Fig.^Jwe display the same information represented as the 
the one electron removal and addition spectral function. 
For non-interacting electrons this represents the occupied 
(left) and unoccupied (right) states. 
In the absence of electron-electron interactions, here rep- 
resented by the on-site Hubbard repulsion parameter U, 
the Fermi energy would be located in the middle of the 
bonding band, resulting in a metallic conductor. In 
Figs.[| and [| we have adopted U = 4 eV. The model now 
predicts a gap of order Eqt — 1 eV. The Fermi energy 
is located within this gap. The fact that these materials 
are insulating therefore is associated with the large on- 
site Hubbard interaction. 

The excitation of an electron across the gap involves a 
change of occupancy of two of the rungs: The final state 
has one empty and one doubly occupied rung. It is im- 
portant in this context, that the two electrons | 3,1 Li? > 
are in a correlated state: In the limit U — * oo one electron 
is located on the lefthand V-atom and the other on the 
righthand V-atom. 



V. DISCUSSION OF THE EXPERIMENTAL 
SPECTRA 

The 0.9 eV peak marks the fundamental gap of the 
optical spectrum. The interpretation of this peak is still 
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FIG. 5. Occupied and unoccupied density of states ob- 
tained within the Heitler-London model for U — 4eV (a), 
and the schematic diagram of the relevant energy levels of 
the V2O clusters forming a rung (b). L and R refer to or- 
bitals associated to the left and right V sites, respectively. 



subject of a scientific controversy. Several interpretations 
have been put forward 

1. Transitions between linear combinations of \L3d_— 
states oLthe .two V-ions forming the rungsuli-ilLil. 
In Refs.cl andll 4 ] even and odd combinations were 
considered. The 0.9 eV peak in a a (uj) (peak A) 
would then correspond to the transition from V-V 
bonding to antibonding_cpmbinations on the same 
rungt il (Fig.|b). In Refill this model was extended 
to allow lop-sided linear combinations of the same 
orbitals, so that the 0.9 eV peak is a transition be- 
tween left- and right-oriented linear combinations. 

2. On-site d-d transitions between the crystal field 
split levels of the V-ionsLH Because V 5+ has no 
occupied 3d-levels, such processes involve the V 4+ 
ions. 



bond- 



3. Transitions between the on- rung V 3d^ y 
ing combination and final states of d xz and d 
charactered. 



yz 



Optical transitions having values below 2 eV were also 
seen in VgOi3 and VO2. In V2O5 they haiie very small 
intensities, and were attributed to defectsHl. The last 
two assignments are motivated by the fact that in a'- 
Nai-zCa^X^Os the optical selection rules allow on-site 
d-d transitions by virtue of the low point-symmetry at 
the V-sites. To determine which one of these assign- 
ments is true, we have measured the doping dependence 
of the 0.9 eV peak in a a (v) in Ca-substituted a'-NaV205 
(Fig.||a). Because Ca is divalent, substituting Na with Ca 
has the effect of increasing the average density of V 4+ 
ions. 

In the local d-d scenario the intensity of the on-site V 
3d-3d transitions would be proportional to the density of 
V 4+ ions. As a result the intensity of on-site d-d tran- 
sitions increases upon substituting Na with Ca. In Fig. 
[| we display the experimentally observed doping depen- 



dence of the intensity together with the theoretical ex- 
pectation within this scenario. Clearly the experimental 
intensity of the 0.9 eV peak in Fig. || behaves opposite to 
the expected behaviour of dd transitions. This definitely 
rules out item number 2 of the above list. _ 
This also rules out item number 3 presented in Ref.tj: 
if the transition from the V 3d+ y bonding combination 
to the d X z and d yz orbitals involves mainly transitions 
among orbitals at the same site, the same argument as 
for item 2 applies. If it involves mainly transitions be- 
tween molecular orbitals formed by different sites on the 
same rung, the transition from the V 3d xy bonding com- 
bination to the antibonding V 3d~ y would still be the 
dominant transition. 

To explain the intensity decrease of about x% upon sub- 
stituting x% of the Na ions with Ca (see Fig.||), let us 
have a look at the many-body nature of the ground state 
and the excited state properties of NaV20s. The rungs 
which are occupied with two electrons due to the Ca- 
doping will be in the many-body ground state (see Table 
|) \ X LR). Due to the dipole selection rules, the optical 
excitation with the electric field along the rung-direction 
is I 1 LB) — ► \ l O). The energy to make this excitation is 
hU + a/<7 2 /4 4- 4t 2 _. Hence the effects of Ca doping are 
(i) to remove the peak at Eqt ~ l&V for the rungs re- 
ceiving the extra electron, and (ii) to place a new peak 
at an energy U ~ 4eV. Hence, the observed x% decrease 
of intensity of the \B) — ► \A) transition peak for the x% 
Ca doped sample is in excellent quantitative agreement 
with the expected value. 

Using the tightbinding f-sum rule J a(uj)duj = 
(ed/h) 2 irt ± /(2V)(d{ a d Ra . + HC) the intensity of the 
\ X LR) — > ^O) peak relative to the 0.9 eV peak of the 
singly occupied rungs is (assuming A = for a doubly 
occupied rung): 
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With the parameters relevant to NaV20s this implies 
that the \ X LR) — > ^O),] 1 !^) transitions have factor 2- 
4 smaller spectal weight than the \B) — > \A) transition. 
Hence we conclude that only the assigment of item num- 
ber 1 is consistent with our data: The leV peak in cr a (uj) 
is the on-rung \ 2 B a >^ \ 2 A a > transition with an exci- 
tation energy Ect = y 1 At\ + A 2 . 

The 1.1 eV peak in <7b(u;) (peak B) involves transitions 
between neighboring rungs along the ladder. In the non- 
interacting model (U = 0) this would correspond to a 
Drude-Lorentz optical conductivity centered at u> = 0, 
with a spectral weight J ab(u})duj = (e//T-) 2 t||7r&(2ac)~ 1 . 
As a result of the correlation gap in the density of states, 
indicated in Fig. ||a, the optically induced transfer of 
electrons between neighboring rungs results in a final 
state with one rung empty, and a neighboring rung dou- 
bly occupied, in other words, an electron hole pair con- 
sisting of a hole in the band below Ej?, and an electron 
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in the empty state above Ep indicated in Fig. ga. This 
corresponds to the process 

2\ 2 B a >^\ 3 - 1 LR>+\0> (5) 

Note that the final state wavefunction is qualitatively 
different from the on-rung bonding-antibonding excita- 
tions considered above, even though the excitations are 
close in energy: it involves one rung with no electron, 
and a neighboring rung with one electron occupying each 
V-atom. The energy of this process is approximately 
2t± + SV, where SV represents the increase in Coulomb 
interaction by bringing two electrons together on the 
same rung. Since the distance betweeen the electrons 
changes from about 5.0A to 3.4A, and taking into account 
a screening factor e « 6, we estimate that SV ~ 0.2eV. 
This value of SV corresponds closely to the difference in 
peakpositions along the a and b directions. According to 
this interpretation the absorption at 1.1 eV along b corre- 
sponds to the creation of a free electron and hole, capable 
of carrying electrical currents. The on-rung excitation at 
0.9 eV along the a-direction is a localized (charge neutral) 
excitation, in other words an exciton. In this case the en- 
ergy of the exciton involves the states of a single electron 
only, whereas the free carrier states involve many-body 
interactions. 

Doping with Ca creates doubly occupied rungs, whose 
ground energy is not 2\ 2 B a > but \ X LR). Consequently, 
the electrons on these rungs will not be involved in the 
processes of Eq. |5|, thus decresing the intensity of the B 
peak upon doping, as seen from Fig. ^b. 
The room temperature crystal structure has four V- 
atoms per unit cell, organized in ladders with up and 
down oriented apical oxygens alternating along the a- 
direction, resulting in a double degeneracy of the elec- 
tronic states discussed above. The coupling between ad- 
jacent ladders lifts the degeneracy of these states, result- 
ing in a "Davidov" splitting of the peaks A and B. This 
can create the two additional "shoulders" in u{lo) at 1.4 
and 1.7 eV fop-peak A and B respectively. 
With ARPESB it has been observed that an energy of 3 
eV separates the V3c? band from the 02p. We therefore 
attribute the peaks at 3.3 eV in a a (u) and the peak at 4 
eV in <Jb{ui) to transitions of the type 

\ 2 B a >^ |Vii? > (6) 

where the 2p hole is located on the oxygen on the same 
rung for peak A, and inbetween the rungs for peak B. 
This is furtherj-supported by previous optical measure- 
ments on V2O5II3 showed a peak at about 3 eV. In V2O5 
all V-ions have a formal V3d° configuration, hence the 3 
eV peak can not be attributed to d-d transitions. How- 
ever, the 02p ~^V3d transitions should appear at approx- 
imately the same photon energy as in NaV20s, which fur- 
ther supports our assigment of the 3 eV peak in NaV2C>5 
to 02p ~^V3d transitions. 



VI. CONCLUSIONS 

In conclusion, we have measured the dielectric function 
along the a and b axes of Ca a; Nai_ a ;V205 for x=0, 0.06, 
0.15 and x=0.20. The 0.9 eV peak in a a (u>) was shown ex- 
perimentally to be a bonding-antibonding transition in- 
side the V2O rung and not a vanadium d-d transition due 
to crystal field splitting. We identified the 3.3eV peak in 
(T a (u>) as the transition from the oxygen orbitals to the 
antibonding one. This stxaacly supports the notion, pre- 
viously expressed in RefsaoO that NaV205 is an insula- 
tor due to a combination of three factors: A crystal field 
splitting, an on-site Hubbard interaction, and an on-rung 
bonding-antibonding splitting of the two Y3d xy orbitals, 
each of which is large compared to the inter-rung hopping 
parameter. 
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Definitions 



| = ^l + [A/2^] 2 + [A/2t x ] £ CT = VI^T^ 

TLi?) = yi(4 T 4i - 4i4 T )l°> * = \/t / 2 /4 + 4^- [7/2 

l 1 ^) = V /|(4 T 4 1 +4 T 4i)|0> if = ^1 + [U/4t ± ] 2 + [t//4t x ] 

TO> = V / |(4 T 4i -4 T 4|)|0> ft = + [*/2A] a + [JT/2A] 

TABLE I. Eigenstates and energies of a rung with N e =0, 1 and 2 electrons. The N e =2 state vectors and energies were 
derived for the limit A <C U. 
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